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NEUTRON THERMAXIZAT ION STUDIES IN HYDROGENOUS MEDIA 
Doilnitsyn Ye.Ya. , Novikov A.G., Ilyasova G.A., 

Fanarin M.V., Smelov V .V • , Stakhanov I.P. » Ste- 
panov A.S., Stupak A. I. 

At present the most important hydrogenous neutron 
moderating materials in reactor technics are water and some 
metal hydrides. 

When neutrons are thermalized in water the effects 
of atomic chemical binding and of the molecular diffusion 
of liquid are revealed. 

When neutrons are moderated in metal hydrides ene- 
rgy exchange between a neutron and moderator goes on by 
portions defined by the oscillation levels of a proton 
which is rigidly bound with metal atoms in a lattice. The 
total cross-section of metal hydrides exhibits large effect 
due to rigid binding of a proton with metal atoms in a lat- 
ticeflj . 

To reveal some peculiarities in forming thermal neu- 
tron Bpectra in these media the experiments on neutron 
thermalization in water and lithium-7 hydride were perfor- 
med and the total cross-sections of lithium-7 hydride, zir- 
conium hydrides with different H-content were measured. 
Experimental results are compared to calculations ta ki ng- 
into account specific features of the moderators studied. 

I. Neutron Thermalization in Water and Boric Acid Solutions 
When studying neutron thermalization in water the 
following problems were considered: 

1) The process of thermal neutron spectrum fo rmin g at dif- 
ferent distances from a source was studied. 

2) Relaxation length of neutron gas temperature was found 
directly from experiment and the effects observed were 
compared with calculations for the given experimental con- 
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ditions. 

3) Absorber concentration influence on thermal neutron spe- 
ctrum forming in infinite medium with uniformly distributed 
sources was studied by the comparison of experimental re- 
sults with calculations performed with a model taking into 
account chemical binding strength [2] ,[3] . 

Experimental data necessary to solve the above men- 
tioned problems were obtained by means of space-energy di- 
stribution measurements of thermal neutrons in water pro- 
duced by spherical source with its spectrum being in one 
case nearly Fermi-distribution and in the second case- 
Maxwellian distribution with T = 760°K. 

§ 1. Experimental Method 

Thermal neutron spectra in water at different dis- 
tances from the source were measured by time-of -flight me- 
thod. 

The arrangement used for "these purposes is shown in 

Pig. I. 

Cylindrical tank with 500 mm in diameter and 600 mm 
high mounted near the graphite— water power reactor shield- 
ing was filled with distilled water. Collimated neutron 
beam passed through tine cylindrical tube ( d = JO mm) onto 
the lead sphere O, = 20 mm) placed in the centre of the 
tank filled with water. This sphere isotropically scatter- 
ing the incident neutrons imitated approximately a spheri- 
cal source. Beam hole placed at 70° to the primary beam 
direction was a tube of rectangular cross-section I5x30mm. 

Thermal neutron beam distribution around the sphere 
measured by indium indicator activation only in the direct 
direction had some assymmetry ( -^20%) which disappeared al- 
ready completely in the region of the beam hole. 

Neutron spectra, set up in water, were studied by 
means of the mechanical selector which was described in [4] . 
The energy resolution in these measurements was ~ 40 psec./m. 
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Beam hole, mechanical chopper and detector being fixed re- 
lative to one another, 'they were able to move along the 
beam hole direction that; permitted to measure spectra in 
water at different distances from the sphere. The data obtai- 
ned by means of indium indicators were used for the norma- 
lization of neutron spectra according to -their amplitudes 
and for the introduction of a correction on the gradient 
effect . 

In all the figures the neutron spectra are represent- 

e ! as! 

where —■ - neutron density in the unit interval of 

Oft: 

time-of — flight ; 
t - time-of— flight in ysec/m. 

Source spectrum, i. e. the spectrum of neutrons 
scattered by the sphere, measured without water in the 
tank, is shown in Fig. 2. 

In water the measurements were carried out at the 
following distances from the sphere: =0; I0mm;20mm; 

50mm; 65mm ( % - the distance from the sphere centre to 
the point investigated) . At each of these distances the 
spectrum was measured twice: without Cd and with Cd 
placed at the beam inlet to the tank. Experimental curves 
taken with Co/ in the primary beam represent neutron spec- 
trum from the source with Fermi distribution, as it is 
seen from Fig. 2. Experimental, curves M ^ for the given dis- 
tance as the difference between two measurements ( without 
(Co! and with Cof in -the primary beam) represent neutron 
spectrum in water caused by the thermal component of the 
reactor beam which ( as it is seen from Fig. 2) is nearly 
Maxwellian with T H = 760°K. 

§ 2. Thermal Neutron Spectra from the Source 
with Fermi Distribution 

Space-energy neutron distribution from the source with 
Fermi spectrum is shown in Fig. 5* In all the points studied 
( including the measurements close to the sphere) thermal 

iu 
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range of the spectra is satisfactorily described by Maxwel- 
lian distributions at temperatures slightly exceeding the 
temperature of the medium and coinciding within experimen- 
tal errors. 

The results of this experiment are compared with 
the calculation of the thermal neutron space-energy distri- 
bution in water carried out for spherical geometry. The cal- 
culation was performed in -approximation in the energy 
range 0 < E< 0.67 ev divided into 15 groups . Thermalixation 
effects were accounted for on the basis of the work by Turohin 
V.F. 5 . At E > 0.67 ev (source region) the neutron spec- 
trum was assumed to be a Permian one. Experimental curve 
obtained with the help of sodium indicators ( E = I.46evj 
was used as spatial distribution of the source. The details 
the method used in this calculation may be found in f 37 
[5] . 

Experimental and calculational curves were normali- 
zed according to distribution maximum IV. 

From Fig. 3 the results of experiment and calc ul ation 
are seen to agree sufficiently well. Calculational curves 
practically coincide with the corresponding Maxwellian di- 
stributions in thermal regions. Therefore both experimental 
and calculational data allow to make a conclusion that with 
Permian spectrum source thermal neutron energy distribution 
in water (in the range £ 4-0. X5ev) does not depend on the dis- 
tance to the source and can be described at any point of 
the medium by the thermal neutron spectrum for infinite 
medium with distributed sources [3] ( accurate within a 
B m a1 ^ correction on the localized neutron leakage, § 4). 

8 3* Thermal Neutron Spectra from the Source with Maxwel- 
lian Distribution Temperature 760 °K. 

Space -energy neutron distribution in water created 
by the source with the spectrum being a thermal component 
of the reactor beam ( difference curves) is shown in Fig. 4. 
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Experimental curves were determined to be close to Maxwel- 
lian distribution at *t- — ^ 20 with temperature =3CW-°K 

( the temperature averaged over the measurements at three 
distances: 20 mm- 303 ± I5°K; 50mm- 306± I5°K; 65 mm- 302 * 
I5°K). Therefore to begin with the distance ~ 20 mm from 
the sphere the spectrum of this source does not influence 
the neutron energy distribution within experimental accur- 


In the range ^-^o<20 mm which may be called a transi- 
tional one the neutron distributions measured are of essen- 
tially different character. These curves were compared with 
the calculation based on the assumption that the process 
of setting up medium spectrum goes through Maxwellian- -like 
forms with the temperature depending on the distance to the 
source. Kinetic equation was solved by Grad method [6] by 
Emit symmetrized polynomial expansion of the neutron di- 
stribution function. Expansion parameter, neutron tempera- 
ture T , was chosen so that the distribution function was 
sufficiently well approximated by the first two members: 

. s P J t > co 

where is the second moment of the distribution function. 

By integration of kinetic equation over the veloci- 
ty space we get the system of equations forT, ft (neutron 
density) and J ( neutron current): 

#£ + f! i-rn^O; „ 

+jrnT= '=r ~ n - 1 ; < 2 > 

9 t 3 BX t L I 

*7 _ k _J_ 9 ftT , 

J L ~ jtc i a • 

where j - me ditun temperature, a - absorption coefficient, 
trL- neutron mass. 

Kinetic coefficients X and T were determined in 
the following way: 

nJ -L l o/ x 6J (x^x)xx' J ■ 


-i- Jolx 67 (x'-+x ) xx '*] , 
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for 2-t 0 = 5 mm, 7"= 360±20°K } for =10 mm, T* =3I7-I5°K 
( Fig. 5). Using the expression (4) and the temperatures 
obtained for different distances it was found that L* = 

= I±0.2 cm . Neutron temperature relaxation length calcula- 
tion on the basis of the model proposed in [2 ] gives for 
our case /, = I.I cm . The dependence of neutron temperatu- 
re on the distance to the sphere calculated according to 
(4) with L =1 cm is shown in Fig. 5 • 

From Fig. 4 it follows that difference curve approxi- 
mation by Maxwellian distributions may be considered t^ be 
satisfactory at the distances where the spectrum becomes 
nearly asymptotic one but this approximation is not valid 
for the region where the source proximity causes the marked 
distribution function anisotropy. A more precise descript- 
ion of the experimental curve shape over the whole transi- 
ent region may be obtained by representing these curves as 
superposition of source spectrum ( Maxwellian distribution 
with T= 750°K) and of asymptotic medium spectrum (Maxwell- 
ian distribution with T =300^K) with the amplitudes depen- 
ding on the coordinates. Similar approach was already repea- 
tedly used in a number of calculational and. experimental 
works, for example, in [?] »[&] • Difference curve expan- 
sion in two such components is shown in Fig. 6. 

Spectrum character at energies ^ O.OI ev was assu- 
med to be completely defined by Maxwellian distribution 
with T =300°K. So asymptotic components proved to be found. 
By subtracting them from difference spectra the curves si- 
milar in shape to Maxwellian distribution with T =750°K 
but with different amplitudes were obtained for all the 
three distances ( 0.3mm and IOmm). 

Densities of the "hot" and "cold” components found as a 
result of this expansion as functions of the distance 
from the source are presented in Fig. 7 • 

An attempt to obtain similar curves by calculation 
was made . Neutron distribution function was represented 
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as the superposition of two components each of which de- 
pended only on one of the temperatures. The calculation 
method used did not differ much from the above mentioned 
one. 

The calculation was performed for two cases: for 
the medium without absorption and with slight absorption. 
From Fig. 7 which represents the comparison of these calcu- 
lations with experiment it is difficult to give preference 
to any of the calculation curves. At great distances, how- 
ever, experimental points for"cold M neutron density sire 
nearer to the calculation data taking into account the 
absorption. 

§ 4. Neutron Spectra in Boric Acid Solutions 

When measuring thermal neutron spectra in boric acid 
solutions the experimental arrangement described above 
was used. Neutron flux from the reactor at the tank inlet 
was always covered with cadmium, and the butt-end of the 
beam tube was placed above the sphere that resulted in 
small gradient of scalar flux in the direction of beam 
direction. 

When comparing the results obtained by the calculat- 
ions for infinite medium with uniformly distributed sour- 
ces it is necessary to introduce a correction for the lo- 
cal diffusion neutron leakage. In the region near the sou- 
rce this correction at first sight seems to be considerable. 
However, the distribution type of neutron flux created by 
our source spectrum is such that beam extraction region 
-I* Icm) is close to the bend point of the curve 
describing the distribution of this flux. Local leakage 
is revealed here as slight additional ab sorption (^r< 0) 
accounting for - 15 % from the total absorption for pure 
water and practically does not affect the experimental 
results with additionail absorber. 

Experimental results [3] and their comparison with 
calculations performed for infinite medium with unif ormly 
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distributed sources [3] are presented in Fig. 8. When plot- 
ting the data^w/7^ as a function of absorber concentration, 
additional absorption caused by local diffusion leakage is 
accounted for. In internal region experimental curves satis- 
factorily agree with calculation spectra which prove to be 
nearly Maxwellian distributions with temperatures shown in 
the same figure. 7^ Rvalues as a function of absorber con- 
centration in water are compared in Fig. 9 with calculations 
performed in 2 3 with accounting for chemical binding for- 

ces (curves I and 2) and according to the gas model (curve3)* 

It mush be pointed out that the results of the last 
calculations described in the paper presented to the present 
conference [2] agree better with our data. 

The whole set of the experimental values 7^/7^, cannot be 
explained on the basis of gas model. 

II. The Study of Neutron Hydride Moderating Properties 

§ I. Total Cross-Sections of Lithium-7 and Zirconium 
Hydrides 

Total cross-sections of metal hydrides were measured 
by means of mechanical neutron selector described in ElOl [ IlJ . 
Spindle-shaped slits in the rotor and the system of syncro- 
nously rotating rotors ensured negligible background in the 
measurements at flight path 130 m with resolution 3 jusec/m in 
thermal region. 

Figure 10 shows the trend of the total cross-section 
curve for the neutrons scattered by hydrogen bound in lithium-7 
hydride according to the data of measurements of lithium-7 
hydride total cross-section from 0.02 ev to 2 ev. 

At neutron energy above 0.2 ev scattering cross-section 
obtained experimentally follows the law f 12 J [ 1 3 j : 

= (5) 

where k y -average proton kinetic energy, p - proton mass 
to neutron mass ratio, ^ - cross-section of scattering 
by a free nucleus. 

Calculational cross-section values for neutron scatter- 
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xng by litium— 7 hydride obt ained in [2] with regard for lev- 
els at 0.1 ev, 0.15 ev and acustic oscillations satisfacto- 
rily agree with, experimental data. In Pig. 10 the dotted 
curve represents the calculated values in the range of 0,1 
ev level where "there is a slight difference as compared to 
the experimental results. Good agreement between experiment- 
al and calculational data is observed in all the rest of the 
range. Zirconium hydride total cross-section was measured 
for zirconium hydrides with different H "content ; 

• The results of the measurements show the ag- 
reement between the total cross-sections of zirconium hyd- 
rides within experimental accuracy (4%) and with the data 
published in [I] . 

No displacement of oscillation levels of H in the 
lattice was observed for zirconium hydrides with different 
H - content. 

Neutron scattering cross-section in the range 
E > 0.5 ev is described by the relation (5). 

§ 2. Neutron Thermalization in Lithium-7 Hydride 
Neutron thermal iz at ion in lithium-7 hydride was 
studied by the method described above in this report (§4). 

The same arrangement, identical geometry and similar meth- 
ods of experimental data analysis were used in these expe- 
riments. An assembly from lithium— 7 hydride had cylindrical 
form with d =4-50 mm and H = 520 mm out of disks with 
d - 47mm and H =5mm . Disks from hydride (canned with al- 
uminium 0.1 mm thick) were packed up so as to reduce the 
escape of neutrons through the assembly slits ( in every 
hydride layer the disks are displaced so as to prevent 
neutron leakage without moderation)* The lead sphere 
( d -50mm) is placed in the centre of the assembly. Neut- 
ron beam from the reactor went through the reentrant hole 
(of square cross-section with side =55mm) onto the lead 
sphere. Cadmium was placed on its way before the assembly. 
Neutron beam went out from the hole butt-end above the 
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sphere centre. Extraction tube was of rectangular cross-sec- 
tion 20x35 mm • 

The results of measurements of thermal neutron densi- 
ty are given in Fig. II a s a function of time-of flight in 
jusQt./m. Statistical accuracy of — determination was 3%. 

To a first approximation thermal neutron spectrum 
in lithium-7 hydride may be described by Maxwellian distri- 
bution with temperature T =33O°+I0°K . The comparison of 
these results with corresponding spectrum measurements in 
water, at the same temperature and with equal absorption 
per hydrogen atom shows the neutron spectrum setting up in 
lithium-7 hydride to be more hard than that in water. To 
evaluate more obviously the hardness of the spectrum obtai- 
ned in lithium-7 hydride, Fig. 12 shows the ratio T H /^ob- 
tained for lithium-7 hydride and the calculational data 
for water [2] . In the same figure calculation curve for 
lithium-7 hydride taken from the same work is presented. 
Neutron absorption for lithium-7 hydride plotted on the 
abscissa axis takes into account the effective cross-sect- 
ion of impurity Li^. From figure 12 good agreement between 
experimental data and calculated data obtained in [ 2 ] la seen. 

The comparison of ^ A>for lithium-7 hydride and wa- 
ter reveals that thermal neutron spectrum setting up in 
hydride is close to the spectrum in water poisoned by ab- 
sorber according to y /V law with absorption in water being 
nearly twice the absorption in lithium hydride. Thermal 
neutron spectrum in lithium-7 hydride Is softer as compared 
to neutron distribution setting up in zirconium hydride at 
the same absorption. 

CONCLUSION 

The experiments on the studies of the moderating 
properties of water, lithium-7 and zirconium hydrides and 
their comparison with calculations show the experimentally 
observed effects to agree well enough with the calculations 
performed with the corresponding differential cross-sectionr 
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Fig. I. Experimental 


arrangement diagram* 


Fig. 2. Source speotmm* 



source with Farai distribution. 
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